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Abstract 

Cold atoms, loaded into an optical lattice with double-well sites, are considered. 
Pseudospin representation for an effective Hamiltonian is derived. The system in 
equilibrium displays two phases, ordered and disordered. The second-order phase 
transition between the phases can be driven either by temperature or by changing 
the system parameters. Collective pseudospin excitations have a gap disappearing 
at the phase-transition point. Dynamics of atoms is studied, when they are loaded 
into the lattice in an initially nonequilibrium state. It is shown that the temporal 
evolution of atoms, contrary to their equilibrium thermodynamics, cannot be de- 
scribed in the mean-field approximation, since it results in a structurally unstable 
dynamical system, but a more accurate description is necessary taking account of 
attenuation effects. 
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I Introduction 



Degenerate cold gases, Bose [1-8] as well as Fermi [9,10], possess many interesting prop- 
erties. Loading cold atoms into optical lattices yields highly controllable systems that 
could be employed for a variety of applications [11-15]. Recently, a novel type of optical 
lattices has been experimentally realized, each site of which is formed by a double-well 
potential [16-21]. 

The phase diagram of cold bosons in double-well optical lattices has been studied 
in Refs. [22,23], where the Hubbard model is used and the main attention is payed to 
the peculiarities of the superfluid-insulator phase transition in such lattices, as compared 
to this transition in the standard single- well optical lattices [11-14]. In Ref. [24], the 
superfluid phase is studied, aiming at finding a type of superfluid with the broken time- 
reversal symmetry. 

The aim of the present paper is to study a different regime that can be realized in 
the double-well optical lattices. We consider the lattices that are in the insulating state, 
far from the boundary of the insulator-superfluid transition. In that state, the jumps 
of atoms between different lattice sites are suppressed, while the tunneling between the 
wells of a single-site double well can be of principal importance. We show that in such 
insulating double-well optical lattices there exists another phase transition, the order- 
disorder phase transition. The investigation of this different physical regime in insulating 
double-well optical lattices, far from the superfluid phase-transition, but close to another, 
order-disorder phase transition, is the main motivation for the present work. 

The order-disorder phase transition, and the related ordered and disordered states, 
can be the most clearly realized for the case of one atom per a double well. Therefore, 
we consider exactly this case. An additional argument for considering the lattices with 
one atom per site is that the optical lattices with small filing factors, such as one or two, 
seem to be good candidates for quantum information processing [11,19]. 

The advantage of dealing with the case of one atom per a double well is twofold. 
The most important, as is emphasized above, is that this is the setup allowing for a 
clear realization of the order-disorder phase transition. Another, technical, convenience 
is that in this case the pseudospin representation can be invoked. Some types of spin 
Hamiltonians using cold atoms in optical lattices can be met in literature. It is straight- 
forward to obtain a spin Hamiltonian for spinor condensates [25]. The Hubbard model 
for a system of two- component bosons can be reduced to a pseudospin representation by 
a second-order perturbation theory in the tunneling parameter [26-28]. In our case of 
an insulating double-well optical lattice, the pseudospin representation can be introduced 
without perturbation theory, by means of an exact canonical transformation, similar to the 
canonical transformations used for deriving the spin representations for superconductors 
and ferromagnets [29,30]. 

The key point of realizing the ordered and disordered states in a double-well optical 
lattice is the existence of sufficiently long-range atomic interactions. Such interactions of 
dypolar type arise between polar molecules [31], Rydberg atoms [32], and between atoms 
with large magnetic moments [33]. The existence of long-range atomic interactions is 
another point making our consideration principally different from the earlier theoretical 
works [22-24] on the double-well lattices. 
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Characterizing collective atomic states in the double- well lattices, we, first, describe 
their equilibrium properties. However, we keep in mind that in experiments atoms need 
to be loaded into a lattice, and their initial state after the loading may be nonequilibrium. 
If so, how then atoms would relax to their equilibrium state? Another way, when atoms 
can happen to be in a nonequilibrium states, is if the lattice parameters are varied after 
the atoms have been loaded. If this variation is sufficiently fast, atoms again occur to 
be in a nonequilibrium state, from which they should relax to an equilibrium one. In 
order that our consideration of the insulating double-well optical lattices would be more 
complete, and keeping in mind the possibility of realizing nonequilibrium states, we study 
not only the equilibrium properties of such lattices, but also the relaxational dynamics of 
atoms from initially nonequilibrium states to their equilibrium. 

The paper is organized as follows. First, we derive an effective Hamiltonian for the 
considered insulating system and show that it allows for a convenient pseudospin repre- 
sentation (Sec. II). Then we study the equilibrium thermodynamics of the model, which 
exhibits the existence of two phases, ordered and disordered (Sec. III). Collective excita- 
tions, corresponding to pseudospin waves, are described in Sec. IV. Dynamics of atoms, 
loaded in the lattice in an initially nonequilibrium state, is considered in Sec. V. Lyapunov 
stability and structural stability of solutions to the equations of motion are investigated 
in Sec. VI, where the principal importance of taking into account attenuation effects is 
demonstrated. The main results are summarized in Sec. VII. 

II Model Hamiltonian 

We start with the general form of the energy Hamiltonian 



where ip(r) = ip(r,t) is a field operator, in which the time dependence, for brevity, is 
omitted. Keeping in mind the case of an insulating lattice, with the unity filling factor, 
the statistics of atoms is not as important. So, atoms can be either bosons or fermions. 
In the lattice Hamiltonian 



the lattice potential V(r + a«) = V(r) is periodic over the lattice {a^ : i — 1,2, ... , Nl} 
and enjoys the double-well structure around each of the lattice sites a,. The interaction 
potential $(— r) = <3>(r), in general, is a sum of the short-range interactions, whose 
strength can be regulated by the Feschbach resonance techniques [2,10,15,34,35], and of a 
long-range interaction of the dipolar type, such that exists between polar molecules [31], 
Rydberg atoms [32], and atoms with large magnetic moments [33]. If, instead of atoms, 
we consider ions, then these exists the long-range Coulomb interaction. The strength of 
the interaction potential is assumed to be such that the intersite interactions, at least 
between the nearest neighbors, cannot be neglected. 

The field operator can be expanded over Wannier functions 



H 






(2) 



c nj w n( r a j) j 



(3) 



nj 
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where n is a band index and j enumerates the lattice sites. Then Hamiltonian (1) trans- 
forms into 



H-V 1 V F mn J r ■ -I- - <T) n i n 2 n 3n4 t „t „ . „ . ( A\ 

ij mn z {j} {„} 

where 

= / <(r - a,) lf £ (r) Wn (r - a,) dr (5) 

and ^jj^j™]™ 4 is the corresponding matrix element of the interaction potential. 
It can be shown (see Appendix A), that 

Eg" = 5 mn E% (6) 

is diagonal with respect to the band indices. Here, 

E% = 5 ij E n + (l-5 ij )J%, (7) 

with 

E n = J w*(r) H L (r) w n (r) dr , J™ = J w* n (r - a^-) H L (r) w n (r) dr . 
Keeping in mind an insulating lattice implies that the intersite hopping is small, such that 

« 1 (< ^ j) • (8) 



Jij 



E n 

Assuming that each lattice site contains just one atom, we impose the unipolarity 
conditions 

^ ] c n jC n j 1 , c n jC n j . (9) 

n 

Taking into account the above consideration, for Hamiltonian (4) we have 



H — E n (^ n jC n j + ^ ^ij C mi C nj C m'jC n 'i , (10) 

nj i^j mnm'n' 

where 

umnm'n' _ ^rrmrn'ra' i ^mnn'm' ("\"\\ 
ij ~ ijji ijij ' V ) 

with the upper or lower sign for bosons or fermions, respectively. 

When a single-well potential is transformed into a double-well potential, then each 
spectrum level of a particle, it would possess in the former potential, splits into two lines 
for the particle in a double-well potential. The spectrum splitting is connected with the 
particle tunneling between the wells of a double-well potential. The splitting magnitude 
depends on the characteristics of the double-well potential and can be regulated in a 
wide range [36]. Thus, for a double- well lattice, one cannot limit oneself by considering 
solely the lowest energy band, but at least two energy levels must be taken into account, 
describing the level splitting and interwell tunneling. In what follows, we take into account 
two lowest energy levels, so that the index n takes two values n — 1, 2. 
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The symmetry properties of the ground-state wave function and of that for the first 
excited state are known [36] to be different. Enumerating the ground state with n — 1 
and the excited state with n = 2, one has 

wi(-r) = Wi(r) , w 2 (-r) = -w 2 (r) . (12) 

Also, both functions w n (r), for n = 1, 2, can be taken to be real. Because of the symmetry 
property (12), the matrix elements of the type as V^ 112 and V? 221 become zero. 
For what follows, it is convenient to introduce the notation 

E = l - (Ei + E 2 ) . (13) 

And let us define the interaction matrix elements 

Mj = \ iylT + V 22 22 + 2V? 21 ) , B iS = \ (V^ + V 2222 - 2V^) , 

C i3 es \ {V 2222 - v£™) , I %3 ee -2V^ 22 . (14) 

The quantity 

n = E 2 - E 1 + Cij (15) 

is the tunneling frequency characterizing the tunneling between the wells of a double-well 
potential. 

The convenience of dealing with the two-level case is that it allows for the introduction 
of the pseudospin representation. The pseudospin operators can be defined as 

Sj = - (AjCij - c 2j c 2 ^j , Sj — - [c[jC 2 j - c 2j cij^j , 



S j = \ ( c l c 2j + 4 c ii) ' ( 16 ) 



which gives 



2. | 1 

c ij c ii = 2 + > c 2j c 2j = 2 ~~ Sj , 



c f ljC2j .s; ; - is] , = s-; + is] . (17) 

The physical meaning of the pseudospin operators (16) can be clarified by introducing 
the left, CjL, and the right, Cj R , location operators 

C jL = ^= (Clj + C 2j ) , C jR EE (cy - C 2j ) , (18) 

characterizing the left or right location of an atom in the left or right well of a double- well 
potential. The pseudospin operators (16), expressed through the location operators (18), 
become 

Sj = - (cj L CjR + Cj R CjLj , Sj= — ~ (cj L CjR — Cj R CjLj , 
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Sj = \ (c] L c jL - c) R c jR ) . (19) 

This representation demonstrates that characterizes the tunneling intensity between 
the left and right wells of a double-well potential centered at the j-site; the operator Sj 
corresponds to the Josephson current between the wells; while S? is the displacement 
operator describing imbalance between the wells. 

Finally, for Hamiltonian(lO), we obtain the pseudospin form 

H = E N + \ J2A i3 - n^S] + £ I^S-SJ _ YshjSlS* . (20) 

The first two terms here are not of the operator type, hence, can be omitted. The third 
term describes the tunneling between the wells of a double-well potential at the j-site, with 
the tunneling frequency Q. The fourth and fifth terms characterize particle interactions, 
with the transverse strength Bij and longitudinal strength i^-. The values of the system 
parameters depend on the properties of particle interactions and on the features of the 
lattice potential. Generally, these parameters can be varied in rather wide ranges. 

In order to illustrate how the tunneling frequency (15) can be varied, we may take the 
double- well potential, in the vicinity of the lattice site &j =0, in the form 



V(r)^V ( r ^) 2 \( r ^ 



+ V H (r y ,r z ) , 



where VH{r y ,r z ) is a harmonic potential in the y- and ^-directions. The tunneling fre- 
quency f2 essentially depends on the parameters of the double-well potential V(r), its 
depth V and the interwell distance r . These parameters enter the variable 



a 



^2mrlV ' 



which varies in the interval < a < oo. For intermediate values of a, the tunneling fre- 
quency can be calculated numerically, while for small and large a, asymptotic expressions 
are available [36] yielding 

n ~ 6V exp (- (a < 1) , 

n ~ V a 2/3 (a > 1) . 

These expressions show that the tunneling frequency Q can be regulated, by changing r 
and Vb, in a wide range < Q < oo. This means that varying the shape of the double- 
well potential, which is achievable in experiments [16-21], one can regulate the tunneling 
properties of the system. 



Ill Equilibrium Phases 

Thermodynamics of the system with Hamiltonian (20), for sufficiently long-range inter- 
actions, can be accurately described in the mean-field approximation [37] corresponding 
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to the equality 

S?Sf = <S«>Sf + S?<Sf>-< 5? >< Sf > (i ^ j) , (21) 

where the angle brackets imply statistical averaging. We introduce the notation for the 
average interaction strengths 



A =~kV Ai " B ^k^ B -" ^k^ 1 "- (22) 

u i¥=J 1+3 i¥=3 



Also, we define the effective tunneling frequency 

Q eff = Q-2B < Sj > (23) 

and the effective mean field 



h eff = y/tf eff + 4/2 < <?J >2 . (24) 
Using the standard methods of dealing with pseudospin Hamiltonians [38], we find 

&eff . ^ I hp 

l eff 



< S ;> = <S| > J-tanhfe), (25) 



L eff 

where T is temperature. The system free energy reads as 

F = NE + y (A - 2B < S? > 2 +21 < S* > 2 ) - NT In ^2cosh j , (26) 

in which N = N L because of the unity filling factor. From the latter equation, one can 
calculate all thermodynamics characteristics. 

We shall be mainly interested in the properties of the average tunneling intensity 

* = W E < ^ > > (27) 

average Josephson current 



^J-£<SJ>> (28) 
L j 



and the average well imbalance 



^J-E<^>- (29) 



The properties of the system depend on the parameters defined in Eq. (22), which 
characterize the effective interaction strength. The parameter A enters the free energy 
(26) additively, hence, it does not play here an important role. The parameter B is 
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composed of the effective interactions Bij given in Eq. (14). The symmetric terms 
V? 222 , and V^ 221 are close to each other, because of which the transverse interactions B^ 
are small, so that the parameter B is substantially smaller than the parameter /. The 
main role among interactions is played by the longitudinal interaction strength /. Here 
we assume that / is positive, which is the often situation for the exchange interactions 
[30,37,38]. Then, as follows from Hamiltonian (20), the nearest spins Sf and S? tend to 
align parallel in order to reduce the system energy. Therefore only the ferromagnetic-type 
order is possible in the system [30,37,38]. 

Generally, the exchange interactions could be negative. Then the longitudinal term 
in Hamiltonian (20) would enter with the sign plus, which would imply that the nearest 
spins prefer to align antiparallel to each other in order to lower the system energy. In 
that case, the antiferromagnetic-type order could be the sole possibility. In the present 
paper, we limit the consideration by a positive parameter I > 0, which means that only 
a ferromagnetic-type order can arise. 

To simplify the formulas, we use the dimensionless quantities, such as the dimensionless 
tunneling frequency 

* - jfs (30 » 

and the dimensionless transverse interaction 

JS JTf (31) 

Also, we define the dimensionless field 

h eff 



'^TTb- < 32) 



which, with the use of the above notations, becomes 



h = yj{u> " bx) 2 + (1 - b) 2 z 2 . (33) 

/,From definitions (27) to (29) and Eqs. (25), we have the tunneling intensity 

uj — bx f h \ .„ „. 

x = — - — tanh — , (34) 



Josephson current in equilibrium 
and the well imbalance 



h \2T 



y = , (35) 



z = z - - - tanh f \ , (36) 
h \2T) 1 K ' 

where temperature T is measured in units of I + B. 

We may notice that Eqs. (34) to (36) are invariant under the replacement 

x — > —x , uj — > — uj , z — > —z . (37) 
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Therefore, without the loss of generality, we can consider only the case with x > 0, u > 0, 
and z > 0. 

Equation (36) shows that there can exist two thermodynamic phases, ordered and 
disordered, when, respectively, 

z 7^ (ordered) , 

z = (disordered) . 
In the ordered phase, the mean well imbalance is nonzero, and one has 

x = u (z^O) , (38) 

with z given by the equation 

i^ ta „h(A) =1 . (39) 

^From Eqs. (33) and (38), it follows that 

h= (l-b)Vuj 2 + z 2 . 

The ordered phase exists, when two conditions are valid, the tunneling frequency is not 
too high, 

< u < 1 , (40) 
and the temperature is lower than the critical temperature 

T c = {1 - b)UJ , (41) 
2 artanh u 

so that 

< T < T c . (42) 

When at least one of conditions (40) or (42) is not valid, that is, either u > 1 or T > T c , 
the ordered phase cannot exist and transfers to the disordered phase. 
In the disordered phase, the mean well imbalance is zero, and one has 

x = tanh {^P^j (z = 0) , (43) 

which defines the tunneling intensity x. 

At zero temperature T = 0, the ordered phase is described by the equations 

x = u , z = Vl - oo 2 (u <1) . (44) 

While the disordered phase is characterized by the expressions 

x = 1 , z = (u > 1) . (45) 

The quantum phase transition between the ordered and disordered phases occurs at the 
tunneling frequency uj — 1. 
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IV Collective Excitations 

Collective excitations in the system represented by the pseudospin Hamiltonian (20), 
correspond to pseudospin waves. The Heisenberg equations of motion for the pseudospin 
operators yield 

fJQV 

-f- = ns? - 2S? J2 hjS* - 2si J2 Bus* , 

^ = -QSf + 2Sf £ B l3 S* . (46) 

For describing collective excitations, we employ the Fourier transformation for the 
pseudospin operators 

ly L k j 

where k pertains to the first Brillouin zone. Similarly, we expand the interaction functions 



% = TT E , B i3 = ±- Y: B h e**** , (48) 



in which a^- = a* — a,. Then, Eqs. (46) acquire the form 

^ = -n*Z + ^ E 5*-^X-p , (49) 

where a k = cr k (t). 

To find the spectrum of collective excitations, we resort to the random-phase approx- 
imation. For this purpose, we look for the solutions to Eqs. (49) in the form 

= < a a k > + Sate-" , (50) 

describing the deviations from the equilibrium values 

<<%>= S k0 N L <S«> . 

Linearizing Eqs. (49) with respect to 5a k yields for the spectrum of pseudospin waves 

e\ = (n - Bx) [n - (B + I k )x] + /(/ + B k )z 2 . (51) 
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In deriving Eq. (51), we have taken into account that 

lim K — I , lim Bi, = B , 

with I and B defined in Eq. (22). 

To consider the spectrum property in the long- wave limit, we use the asymptotic 
equalities 

h * I - - \ £ • a 4J ) 2 , B k ~ - I £ %( k ' ««) 2 , (52) 

i i 

where fc — > 0. We define the gap A by the equation 

A ^ 2 



J + B 

And let us introduce a matrix G^, given by the equality 



uo-bx)(uo-x) + (l-b)z 2 . (53) 



= ( W - - (1 - b)z 2 B tJ . (54) 

Then the long-wave limit of spectrum in Eq. (51) gives 

4^ A 2 + \ EG#^ 3 ) 2 . (55) 

The spectrum is quadratic, which is typical of spin waves. 

Specifying the gap A and the matrix Gij for the ordered and disordered phases, we 
have for the ordered phase, when z ^ 0, 

7^ = z, ^ = (1 - b) (u% - z*B l3 ) , 

and for the disordered phase, with z — 0, 

= y/(u-bx)(u-x) , = (w - ftx)^- . 

The gap disappears at the critical point of the order-disorder phase transition. Approach- 
ing this point from the side of the ordered phase, one has z — » 0, while approaching the 
point from the side of the disordered phase, one has x — > uo. In both these cases, A — > 
at the critical point. 



V Nonequilibrium Loading 

Let us suppose that atoms are loaded into a double-well lattice in an initially nonequlib- 
rium state. For describing their equilibration process, one needs to study the temporal 
behavior of the variables defined in Eqs. (27) to (29) and characterizing the tunneling 
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intensity x = x(t), Josephson current y = y(t), and the well imbalance z = z(t). The evo- 
lution equations for these quantities can be obtained by averaging the operator equations 
(46). For this purpose, one often employs the mean-field approximation. This approx- 
imation does not take into account the relaxation processes due to atomic interactions. 
However, such processes could be important in order that equilibrium would be possible. 
Therefore, we shall resort to a more accurate approximation, called the local-field approx- 
imation. This approximation, suggested by Wangness [39], treats atomic interactions as 
occurring in a local field formed by other particles, which results in the appearance of 
relaxation characterized by the attenuation parameters, 71 and 72 whose values can be 
calculated through the given atomic interactions. The local-field approximation assumes 
the existence of local equilibrium [40], so that at each moment of time the variables tend to 
relax to the locally-equilibrium state. Averaging Eqs. (46) in the frame of the local-field 
approximation [39] yields the system of equations 

dx dy 

— = (1 - b)yz - 72(2: - x t ) , = (w - x)z - 72(1/ - y t ) , 

dz 

— = (bx - u)y - 71(2 - z t ) , (56) 

where the time variable is measured in units of 1/ (I+B), while the attenuation parameters 
71 and 72, in units of I + B, and the local fields are 

uj — bx ( h \ 

x t = — ; tanh — - , y t = , 

h \2TJ ' m 

1-b ( h \ 

« = — *^{2f) ■ (57) 

Here T is the temperature, in units of I+B, corresponding to a would be equilibrium state, 
the type of the latter being defined by the system parameters. The explicit expression 
for h is given in Eq. (33). The evolution equations (56) are supplemented with the initial 
conditions 

x(0)=x o , 2/(0) = 2/0 , z(0)=z , (58) 

which, generally, represent a nonequilibrium state. 

Keeping in mind applications to ultracold atoms, we may set temperature to zero. 
Then, Eqs. (56) possess two types of stationary solutions. One type is defined by the 
fixed point 

x\ = u, yl = 0, z\ = VT^ZP, (59) 
which corresponds to the ordered phase. While another type is given by the fixed point 

x* 2 = l, yl = 0, z* 2 = 0, (60) 

corresponding to the disordered phase. The question remains whether an atomic system 
loaded into the double-well lattice, in an initially nonequilibrium state, would relax to 
one of the stationary solutions. 
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VI Stability Analysis 



The stability of stationary solutions can be studied by the Lyapunov stability analysis. 
To this end, we calculate, in the standard way, the Jacobian matrix associated with the 
dynamical system (56). Then we find the eigenvalues of the Jacobian matrix, denoted 
by A, which give the characteristic exponents. The real parts Re A define the Lyapunov 
exponents describing the stability properties. For simplicity, we set 71 = 72 = 7- 

Accomplishing this procedure for the fixed point (59), we find that it is stable when 
uj < 1. One of the characteristic exponents is exactly 

Ai = -7 • (61) 

The expressions for two other characteristic exponents are rather cumbersome, because 
of which we write them down to the first order in 7, resulting in 

2 — b — uj 2 

A 2,3 ^ - 2(1-6) 7 ± iiUeff ' ( 62 ) 

with the effective frequency 



uj eff = y/(l -6)(1- uj 2 ) (w<l). (63) 

This shows that the fixed point (59) is a stable focus, when uj < 1 and 7 is small. Recall 
that by its definition (31), one always has b < 1. 

The fixed point (60) is stable for u > 1. Then, one of the characteristic exponents is 
exactly the same as Eq. (61). Two other characteristic exponents, to the first order in 7, 
are 

2w - 1-6 . , 

' 3 ~ ~ 2 {u-b) 7± ^//' (64) 

with the effective frequency 



uj. 



eff 



y/(u-l)(u-b) (uj > 1) . (65) 



Thus, the fixed point (60) is also a stable focus for uj > 1 and small 7. 

The value of the tunneling frequency uj = 1 is a bifurcation point related to the 
dynamical phase transition. At this point, the dynamical system is neutral, with the 
characteristic exponents 

Ai = -7 , A 2 = , A 3 = -7 (u = 1) 

for any 7. 

An important observation, following from the above analysis, is that the dynamical 
system is structurally stable only in the presence of the damping parameter 7 > 0. Setting 
the latter to zero leads to the characteristic exponents for the ordered fixed point (59) 

Ai = , A 2 , 3 = ±iw e // (7 = 0, uj < 1) 
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and for the disordered fixed point (60) 

Al = , A 2 ,3 = ±^e// (7 = 0, UJ > 1) 

where u e ff and uj' e ^ are given by Eqs. (63) and (65). These characteristic exponents 
demonstrate that the dynamical system is structurally unstable. 

But the case of no attenuation corresponds to the usual mean-field approximation. 
Hence the above consideration shows that the mean-field approximation results in a struc- 
turally unstable dynamical system, thus, this approximation cannot correctly describe the 
evolutional processes of atoms in a double-well lattice. 

In order to illustrate in an explicit way the difference between the temporal behavior of 
solutions to the evolution equations (56) and the related mean-field approximations, with 
no attenuation, we solve Eqs. (56) numerically. The initial conditions (58) are chosen so 
that to correspond to an initially nonequilibrium state. 

In Fig. 1, the relaxation to the ordered stationary solution (59) is shown for 7 = 1 
and uj < 1. Taking the attenuation parameter 7 equal to one in dimensionless units 
implies that it is of order / + B in dimensional units. For these parameters, the fixed 
point (59) becomes a stable node. Contrary to this, the mean-field approximation, with 
7 = 0, exhibits permanently oscillating solutions. The curves for x(t) and z(t) always 
oscillate around the given initial conditions, which is not a correct behavior. In turn, the 
latter is caused by the structural instability of the dynamical system in the mean-field 
approximation. The relaxation to the disordered stationary solution (60) is demonstrated 
in Fig. 2 for uj > 1 and 7 = 1, when Eq. (60) corresponds to a stable node. But the mean- 
field approximation, with 7 = 0, again exhibits an incorrect oscillatory behavior depending 
on the choice of the initial conditions (58). Again, it is the structural instability that is 
responsible for the incorrect dynamics in the mean-field approximation. 

Full equations (56) can also exhibit the oscillatory behavior of their solutions. This 
happens when the fixed points are the stable focuses, which occurs for 7 <C u>, as follows 
from Eqs. (62) to (65). In particular, if uj 3> 1, than, according to Eq. (65), the oscillation 
period approximately is 2tt/u. This situation is illustrated in Fig. 3. 

VII Discussion 

We have considered a system of atoms in a double-well optical lattice. The case of 
an insulating lattice with the unity filling factor is studied. This type of double-well 
optical lattices is of special interest, being a convenient setup for realizing the ordered 
and disordered states of atoms in a double-well lattice. 

The principal difference of the present paper, as compared to the earlier theoreti- 
cal works on the double-well lattices [22-24,41], where the superfluidity-insulator phase 
transition is studied for atoms with local interactions, is that we consider an insulating 
double-well lattice, with atoms possessing long-range interactions. An effective Hamilto- 
nian in the pseudospin representation is derived. The system exhibits two thermodynamic 
phases, ordered and disordered. The phase transition between the phases can be driven 
either by temperature or by varying the system parameters. For instance, one can vary 
the strength of atomic interactions or the shape of the double-well potential. 
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The spectrum of collective excitations has the form typical of that for pseudospin 
waves. In the long-wave limit, the spectrum is quadratic, with a gap. The latter disappears 
at the phase transition point. 

Nonequilibrium properties of the system are investigated. Physically, the situation 
can correspond to atoms loaded into a lattice in an initially nonequilibrium state. Such a 
nonequilibrium state can also be prepared by disturbing the lattice by external fields. 

The relaxation of the system, from an initially nonequilibrium state to equilibrium, 
cannot be described by the mean-field approximation. This is rather clear from the 
physical point of view, since the mean-field approximation does not take into account 
attenuation effects, hence, is not able to describe the equilibration process in principle. 
iFrom the mathematical point of view, as we show, the mean-field approximation results in 
a structurally unstable dynamical system. To correctly describe the process of relaxation, 
the damping effects, caused by atomic interactions, must be included. This can be done, 
e.g., by employing the local-field approximation. 

The long-range interactions that would be sufficient for realizing the effects described 
above could be of dipolar type, such that occur for polar molecules [31,42], Rydberg atoms 
[32], and quantum gases with dipolar atomic interactions (see review articles [33,43]). A 
good candidate would be the gas of cold atoms 52 Cr, possessing large magnetic moments 
of 6/iB and, as a result, sufficiently strong dipolar interactions [44,45]. Because of the 
existence of several species exhibiting long-range dipolar intreations, the latter can be of 
different strength. In addition, these interactions can be effectively modulated by external 
magnetic or electric fields, thus, tuning the interaction strength in a wide range [46-49]. 

In the dynamics of atoms, relaxing from a nonequilibrium state to their equilibrium, 
an important role is played by the attenuation parameter 7. The calculation of the latter, 
for systems with dipolar interactions, is a well known procedure, described in detail in 
literature on magnetic resonance [50-55]. In the case of particles with magnetic moment 
p, interacting through magnetic dipolar forces, an accurate value of the attenuation pa- 
rameter is given by the expression 7 = pzop 2 /h, where p is the mean particle density 
and Zq is the number of nearest neighbors. This value of 7 is of the order of the interac- 
tion strength. That is why setting, in dimensionless units, 7 to one, as we have done in 
numerical calculations, is absolutely natural. 

Taking, for illustration, the magnetic moment of 52 Cr, equal to p — 6/ib, with pb — 
0.927 x 10~ 20 erg/G being the Bohr magneton, the typical atomic densities in a trap 
p ~ (10 12 — 10 15 ) cm -3 , and accepting for the number of nearest neighbors Zq ~ 10, 
we have 7 ~ (1 — 10 3 ) s -1 . This means that the relaxation time T rd = I/7 is of order 
T re i ~ (10 -3 — 1) s. The lifetime of trapped atomic systems can vary in the range of 
texp ~ (1 — 100) s. Therefore, the relaxation time T re \ in many cases is much smaller than 
t exp , which implies that taking into account the attenuation is of critical importance. The 
simple mean-filed approximation without taking account of the relaxation effects would 
be qualitatively incorrect. 

In conclusion, it is worth noting that the shape of the double-well potential can be 
easily regulated. Hence, it is possible to organize any required sequence of transitions 
between the ordered and disordered states in the double-well lattice during the system 
lifetime. 
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Appendix A 

The eigenproblem for the lattice Hamiltonian (2), 

H L (r)ip nk (r) = E nk ip nk (r) , 

defines the Bloch functions y? n fc( r ) an d the Bloch spectrum E nk . This eigenproblem, 
invoking the relation between the Bloch and Wannier functions 

V?nfc(r) = ]Tu> n (r - a,) e ika ' , 

can be rewritten as 

H L (r) w n (r - SLj) = ^ E nk e ik ^' w n (r - an) , 

ik 

where a^- = aj — a,-. Using this in Eq. (5) gives Eq. (6), with 

ly L k 

Representing the latter in form (7) yields 

k 

and 

L k 



These properties are used in deriving Hamiltonian (10). 
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Figure Captions 



Fig. 1. Dimensionless variables describing the tunneling intensity x(t), Josephson 
current y(t), and well population imbalance z(t) as functions of dimensionless time for 
id — 0.1 and b = 0.5. Initial conditions are Xq = 0.66, yo — 0.75, and z = 0. The case of 
no attenuation (7 = 0) is shown by the dashed curve. The case with attenuation (7 = 1) 
is represented by the solid line. 

Fig. 2. Dimensionless variables x(t), y(t), and z(t) as functions of dimensionless time 
for uj = 1.5 and b = 0.5. Initial conditions are x = 0.33, y$ = 0.5, and z = 0.8. The 
attenuation parameters are: 7 = (dashed line) and 7 = 1 (solid line). 

Fig. 3. Population imbalance for large tunneling 00 = 100, with 7 = 1 and b = 0.5 as 
a function of dimensionless time. Initial conditions are x = 0, y = 0, and zq = 1. 
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Figure 1: Dimensionless variables describing the tunneling intensity x(t), Josephson cur- 
rent y(t), and well population imbalance z(t) as functions of dimensionless time for uj — 0.1 
and b = 0.5. Initial conditions are xq = 0.66, yo = 0.75, and z$ = 0. The case of no at- 
tenuation (7 = 0) is shown by the dashed curve. The case with attenuation (7 = 1) is 
represented by the solid line. 
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Figure 2: Dimensionless variables x(t), y(t), and z(t) as functions of dimensionless time 
for uj = 1.5 and b = 0.5. Initial conditions are x = 0.33, y = 0.5, and z = 0.8. The 
attenuation parameters are: 7 = (dashed line) and 7 = 1 (solid line). 
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Figure 3: Population imbalance for large tunneling u = 100, with 7 = 1 and b = 0.5 as a 
function of dimensionless time. Initial conditions are xq = 0, yo = 0, and z$ = 1. 
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